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Abstract 

We review our proof that in a scaling limit, the time evolution of a quan- 
tum particle in a static random environment leads to a diffusion equation. In 
particular, we discuss the role of Feynman graph expansions and of renor- 
malization. 



1 Introduction 

The emergence of irreversibility from reversible dynamics in large systems has 
been one of the fundamental questions in science since the days of Maxwell and 
Boltzmann. The famous debate about the statistical character of the second law 
of thermodynamics and the related controversy about Boltzmann's Stofizahlansatz 
in the derivation of his transport equation has been very fruitful for physics and 
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mathematics. After Lanford's rigorous justification of the Boltzmann equation for 
a classical many-particle system at short kinetic time scalesIH, the mathemati- 
cal justification of the Boltzmann equation at longer timescales has remained a 
challenge up to the present time. The analogous statement for quantum systems 
remains open even at the short kinetic timescale. 

A related important question is to understand how Brownian motion emerges as 
an effective law from time-reversal-invariant microscopic physical laws, as given 
by a Hamiltonian system or the Schrodinger equation. Kesten-PapanicolaouO 
proved that the velocity distribution of a classical particle moving in an environ- 
ment consisting of random scatterers (i.e., Lorenz gas with random scatterers) con- 
verges to a Brownian motion in a weak coupling limit in dimensions d > 3. In this 
model the bath of light particles whose fluctuations lead to the Brownian motion of 
the observed particle is replaced with random static impurities. A similar result was 
obtained in d = 2 dimensions JH. Recently 0, the same evolution was controlled 
on a longer time scale and the position process was proven to converge to Brownian 
motion as well. Bunimovich and Sinai[5 ] proved the convergence of the periodic 
Lorenz gas with a hard core interaction to a Brownian motion. In this model the 
only source of randomness is the distribution of the initial condition. Finally, Diirr, 
Goldstein and Lebowitz[6 ] proved that the velocity process of a heavy particle in a 
light ideal gas, which is a model with a dynamical environment, converges to the 
Ornstein-Uhlenbeck process. 

Although Brownian motion was discovered and first studied theoretically in 
the context of classical dynamics, it also describes the motion of a quantum parti- 
cle in a random environment, on a timescale that is long compared to the standard 
kinetic timescale (7] |9l. In the following we describe this result and the strat- 
egy of the proof in a bit more detail. Besides the motivation discussed above, the 
random Schrodinger operator that we study is also the standard model for trans- 
port of electrons in metals with impurities, which plays a central role in the the- 
ory of the metal-insulator transition tfTOl ITTi . The outstanding open mathematical 
question in this area is the proof of the extended states conjecture, stating that in 
dimensions d > 3, at weak disorder, the spectrum of such Hamiltonians is ab- 
solutely continuous. Despite much effort, this conjecture has up to now only been 
proven |[T2l[T3l[T4l on the Bethe lattice, which can be interpreted as the case d = oo. 
In a system with a magnetic field, the existence of dynamical derealization at cer- 
tain energies near the Landau levels has been proven recently! 15 ]. 
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2 The problem and the main result 



We consider random Schrodinger operators, both on a lattice and in the continuum, 
in d > 3 dimensions. In this presentation, we focus on the case d = 3. The time 
evolution of the Anderson Model (AM) is generated by 

i^V(0 = Hil)(t), V(0) = with H = --A + XV U on £ 2 (Z d ) (1) 

where —A is the standard discrete Laplacian and the potential is given by V{x) = 
J2 aG %d V a (x), with V a (x) = v a 5 Xta , and v a independent identically distributed 
(i.i.d.) random variables. We assume that m& = E (y™) satisfies 

\fi < 2d : rrii < oo, mi = = m§ = 0, m-2 = 1. (2) 

The continuum analogue of this model is the Quantum Lorentz Model (QLM), 
where H = -±A + AK, on L 2 (R d ), with A the standard Laplacian, V u (x) = 
J Rd B(x — y)dfi UJ (y), where B is a fixed spherically symmetric Schwarz function 
with € supp B, ^ is a Poisson point process on M. d with homogeneous unit 
density and i.i.d. random masses: 

oo 

/ J u = ^ tt r( w )*»»' ( 3 ) 

7=1 

{y 7 (u;)} is Poisson, independent of the weights {u 7 (u>)}. Again, := E„ is 
assumed to satisfy ©. 

Suppose the initial state is localized, i.e. ipQ is smooth. How does the solu- 
tion i/j(t) = e~ ltH i[)o behave for large t ? If A = 0, the time evolution is easily 
calculated in Fourier space: ijj(t, k) = e ~ lte ( k ) Tp {k) , with e(k) = k 2 /2 (QLM) or 
e(k) = X)i=i(l — cos hi) (AM). It is equally easy to see that the motion is ballistic, 
i.e. 

(x 2 ) t = m),x 2 ip(t)) ~t 2 . (4) 

If A 0, one expects either localization, (X 2 ) t = 0(1) for all t, or diffusive 
behaviour (extended states), (X 2 )t = 0{t), depending on A and tfjQ. The localized 
behaviour corresponds to dense pure point spectrum at almost every energy; this 
was proven for large disorder |[T6l[T7l and away from the spectrum of the Laplacian. 
Extended states correspond to absolutely continuous spectrum. As mentioned, the 
latter has been proven lTT2l [T3l [T4l on the Cayley tree for small A > 0. At this time 
there is no proof of existence of extended states in d = 3. For a simpler case, 
namely that of randomness with a decaying envelopping function, i.e. V w {x) = 
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uo x h{x), u! x i.i.d., h fixed, there is a proof |[T8l [T9l that for t] > i and /i(cc) 
as |x| —> oo, H = —A + has absolutely continuous spectrum. 
Our result is formulated in terms of the Wigner function 



W4x,v) = j dy^y^{x + V -)^(x- y -) (5) 

which can be thought of as an analogue of a phase space density (but can be- 
come negative). Its marginals are f W^(x,v)dx = \ip(v)\ 2 and J W^(x,v)dv = 
\i^>(x)\ 2 . Moreover, 



W^,v) = J dx e-^W^v) =^{v -e/2)^(« + e/2). 



(6) 



On the lattice, one has to modify the definition of the Wigner transform slightly 0. 
The kinetic scaling is given by 

■q = A 2 , T = rjt, X = r/x, (7) 

i.e. the microscopic time and space variables both become of order A~ 2 , so that 
velocities remain unsealed. 

Theorem 2.1 

EW J ( r,ri)(*.V) ^ F (X,V,T), (8) 
F the solution of the linear Boltzmann equation 

Ajr(Af,V,r) + (Ve)(V) • V X F(X, V,T) 

= 2nJ dU S(e(U) - e(V)) B{U-V) 2 [F(X, U, T) — F(X, V, T)] (9) 

This theorem was first proven for the continuum for small time T (20], then for 
arbitrary time(2T], and later extended to the lattice case[22|. 
The diffusive scaling is defined by 

£ = X 2+k/2^ x = ex, T = e\ R ' 2 t = \ K+2 t. (10) 

This is long compared to the kinetic timescale: the kinetic variables X and T 
diverge as A — > when X and T are kept fixed, 

X = \- R I 2 X, T = \- K T. (11) 

A first hint at diffusion is that under this scaling X 2 /T = X 2 /T is independent of 
A. The result for the Anderson model is 
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Theorem 2.2 Let d = 3, t/jo £ l 2 (l?) and ip(t) be the solution to the random 
Schrodinger equation with initial condition ipQ. If k > is small enough and 
e = \ 2+K / 2 t then in the limit A — > 0, EW^A- 2 -«y) converges weakly to the 
solution f of a heat equation. 

More precisely: denote (F) E = <S>(E)~ l f dv F(v)6(E-e(v)), where &(E) = 
J dv 8{E — e{v)). Let E £ (0, 3) and Dij(E) = 27r <s>(E) * e ^j e )E, and let f 
be the solution of the heat equation 

^Lf(T,X,E) = Vx-D(E)V X f(T,X,E) (12) 
f(0,X,E) = S(X) {\M 2 ) E . (13) 

Let 0(x, v) be a Schwartz function on R d x M d /27rZ d . Then 

lim V Idv 0(X, v) EWl n - K - 2T) (X, v) 

= / dX fdvO(X,v) f(T,X,e(v)). (14) 
The limit is uniform on [0, To] for any Tq > 0. 

We discuss some of the ideas in the proof of this theorem in Section [3] 

If ipo € C 1 and A is small enough, we have the more detailed error estimate 

Jdvjdt 0(£, v) (£, v) (15) 

= JdtJ $(E)dEe-^ D W^(6(t,-)) E (W^,-)) E +o(l). 

The Boltzmann equation also gives the same diffusion equation in the long time 
limit, but it was itself derived from the quantum mechanical time evolution only 
for shorter timescales. The main difficulty in the proof is to deal with contributions 
that vanish for A — > under kinetic scaling, but that become important under 
the above-defined diffusive scaling. More technically speaking, in the Feynman 
expansion done to analyze the time evolution, most of these terms would even 
diverge under diffusive scaling if we did not renormalize the propagator. 

The allowed values of k are in an interval [0, Kq), where kq is a universal con- 
stant. For technical reasons, has to be chosen very small in the proof. Heuris- 
tically, i.e. ignoring many of the technical complications and assuming optimal 
bounds, one would expect the remainder of the renormalized Feynman graph ex- 
pansion to vanish up to kq = 2, and to diverge for kq > 2. 
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The diffusive scaling leads to a diffusion on the energy shells. A diffusion 
mixing energy shells is expected to start at t = A~ 4 . 

An intuitive way of interpreting the expansion described below is as an ex- 
pansion in the number N of collisions of the particle with the randomly (but stati- 
cally) arranged obstacles represented by the potential. As compared to the previous 
results ll2Tll22l . the main new feature here is that under diffusive scaling, the effec- 
tive number of collisions of the particle diverges. That is, not only is it necessary to 
expand to an order N that diverges as A — > 0, but also the main contribution does 
not come from terms with a finite number of collisions. 



3 Collision histories, Feynman graphs, and ladders 

We discuss some of the ideas of the proof for the example of the Anderson model, 
i.e. the lattice situation. For the detailed bounds of Feynman graphs, the lattice 
leads to a number of complications (H, but for the presentation it is easier. 



3.1 Collision histories 

Let us start with a formal time-ordered expansion, setting Ho = — and ex- 
panding in XV. Then ip(t) = e~ itH t/j = Y^ n >o ^ (n) 0) with 

^W(i) = (-i\) n J dfi n+1 {s)e-' 1SnHo Ve- iSnHo ...Ve-' lsoHo ^ (16) 

where s = (so, ■ ■ ■ , s n ) and 



(17) 



d/u n+ i(s) = J ds . . . ds n S it - , 

[0,00)™ + ! V j = 

Because V = J2aez d Vo. it is natural to split each tp^ further, ip^ (t) = J2a n (*)• 
Every sequence of obstacle labels a n = (ax,... , a n ) € (Z rf ) n represents a colli- 
sion history, and for k 6 {1, . . . , n — 1}, the time variables s& in (fTTT ) are the time 
differences between two subsequent collisions. The delta function in (fTTT ) enforces 
the constraint that these time differences, together with the propagation times sq 
before the first collision and s n after the last one, add up to the total time t. We 
shall discuss convergence questions about this expansion later. 

Our detailed analysis takes place in momentum space, where each V acts as a 
convolution operator, so that 

„ „n— 1 n 

Mt,p n ) = (-iTj dv n+1 (s) / II Wy 1 e~ iSje(pj) Y[ V(pj-Pj-i)$o(pQ). (18) 
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Very schematically, one can represent this as follows, where each of the dashed 
lines represents a factor XV and each of the full lines gets a phase factor e~ ls i e M . 



As the convolution formula shows, one can define a momentum flow in this graph, 
where the momentum change pj — Pj+i flows away through the dashed line. Be- 
fore the disorder average, there is no translation invariance in the system, so every 
scattering at an obstacle changes the momentum of the particle. 

3.2 Disorder average and graphs 

Recalling (©, we have 



E 



n,n' a„ 



4>tJ(t,v-Z/2)^'\t,v + t/2) 



(19) 



Note that there are now two, a priori independent, collision histories, one for ift and 
one for rp. It will be part of the proof to show that, in the scaling limit we consider, 
the only contributions after self-energy renormalization come from the so-called 
ladder graphs, where the two collision histories are identical: n = n' and a n = af n . 

Because the disorder is i.i.d., translation invariance holds for the average, which 
means that momentum conservation also holds for the dashed lines, which for the 
Anderson model simply correspond to a factor A 2 , since the second moment of the 
disorder was normalized to 1 in (0. 

The result can be represented as a graph built of two particle lines, particle- 
disorder vertices, which are joined by disorder lines, and, if the randomness is 
non-Gaussian, disorder-disorder vertices, which correspond to the higher moments 
of the disorder distribution. An example is 



Particle lines get propagators e~ is ' e ^\ interaction lines give factors A 2 , and the 
disorder-disorder vertex of degree four corresponds to a factor m&X 4 . 

It is clear that in the way the expansion was introduced above, one really needs 
the assumption that arbitrary moments, not just the first 2d ones, exist. The ex- 
pansion employed in the true proof contains a stopping rule which avoids high 
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moments, but we shall not discuss this here in more detail. In fact, we shall in the 
following assume for simplicity that the disorder is Gaussian, so that there are no 
vertices of higher degree for the dashed lines, and the average just corresponds to 
a pairing of interaction lines. 

An example of a pairing is as follows 




Note that here, there is a crossing of the two pairing lines in the graphical repre- 
sentation, but there are no vertices in which more than one interaction line enters. 

A special class of pairings are the up-down pairings, where n = n' and the 
pairing corresponds to a permutation a 6 S n : 




The most important term turns out to be the ladder graph, corresponding to 
a = id: 




3.3 Graph bounds 

In the following, we give a brief discussion of bounds of the contributions of indi- 
vidual graphs, restricting to up-down pairings. If one takes a bound in the repre- 
sentation (fT8T >. each phase factor is replaced by 1. This leads to a bound of order 
(Xt) n /nl (where the n! comes from the time ordering implied by the delta func- 
tion in (flTT)). which does not even allow to consider the kinetic scaling where X 2 t 
is fixed. For this reason, the following propagator representation is useful. Let 
7] > 0. Then, inserting the Fourier representation of the delta function, 

„ / n+l \ n+l 

/ d n+1 s 5 [t - V s, ■ TT e~ is ^ 

w +1 ^ k )fi 
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r ( n+l \ n+l 

/. n+l 
da g -ita / d™ +1 S TT e~ iSj<a-e< - Pj - )+iT/ - ) 

2 " i[0,oo)«+l 

/n+l -. 
^ e" iat TT rT — - . (20) 

It is convenient to choose rj = t~ 1 . 

The contribution of a permutation a € <S n , corresponding to an up-down pair- 
ing graph IV, to (6,W$ is 

y G /(r CT ) = A 2n e 2t " y e^-")* 

/ d ^/n0/n& 6(e,ft»X(e,Po) 



j=0 " fe=0 

1 



uj{ qj - f ) - if? « - w(pj - f ) - if? 



n 



n 5 (pj - Pj-i - (9«t(j) - ?<t(j)-i)) • (2i) 

At the moment, w(p) = e(p) € R; later, w will change under renormalization and 
become complex. 

A simple Schwarz inequality separating the dependence on the p, L and that on 
the qi implies that for all a 

\Val(T a )\<Val(T id ). (22) 
The ladder is easy to calculate at £ = 0, and a ladder of length n is of order 

A crucial observation is that the values of graphs with crossings get inverse 
powers of t, as compared to the ladder. This follows from the bound 



/ 



dp-. t~t — — 1 < C| log 77| 3 ,„ I b (23) 

\a - u{p) + irj\ |/3 - uj(±p + q) - ir?| Nil + ?7 



(6 = for the continuum; 1/2 < b < 3/4 on the lattice). |||p||| = \p\ in the 
continuum, |||p||| = min{|p — v\ : V{ € {0, ±7r}} on the lattice. Again, here 
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co(p) = e(p). This motivates why the ladder graph gives the dominant contribution 
under kinetic scaling. However, the number of graphs goes like nl, which cancels 
the 1/n!, hence expanding to infinite order one gets majorants by geometric series, 
which converge only on very short kinetic timescales T. This is the reason for 
the restriction to small kinetic timescales in the first proof [20] of the Boltzmann 
equation for the QLM. 

3.4 Expansions to finite order and remainder terms 

Major progress! 21 1 came from the realization that one can do an expansion to finite 
order with an efficient remainder estimate. A natural way to generate a finite-order 
expansion is the Duhamel formula 

</,(t) = e-' ltH 4> = e- itH °ip + f ds e ^ s)H \Ve-' lsH ^ Q . (24) 

Jo 

Iteration gives 

JV-l 

m = Y,^ n) (t)+^N(t), (25) 

n=0 

where ^ 

x& N (t) = (_i) f dse-^-^W^-^is) (26) 
Jo 

and 

^W(t) = (-iA) n f dfi n+1 (s)e-' ls " H °V . . . V e~ isoHo i; . (27) 



An alternative way of looking at this is via its relation to the resolvent formula 

R z = Rf) + R Z \VR^ (28) 

where R z = (z — H) -1 and R^ = [z — Ho)' 1 . Iteration of the resolvent equation 
and using the Fourier transform gives the above propagator representation directly. 
The Duhamel formula is obtained via 



c 



-UH 



e t V / da e -iat R (29) 



The second crucial ingredient is that one can use the unitarity of the full time 
evolution to reduce all terms to ones where no H appears in the time evolution 
any more: 



*iv(t)|| < / ds 



o 







(30) 
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Thus 



ll*jv(*)ll 2 <*l<M 



ds 



(31) 



The remaining integral over s effectively gives a factor t, which is the price to 
pay for this unitarity bound. To control this factor, one needs exhibit more factors 
i -1 in graphs with several independent crossings, and treat graphs with only one 
crossing explicitly (in the resolvent iteration, the unitarity bound would be replaced 
by ll-Ra+i,,!! < rj~ l ). 

By a Schwarz inequality, one can see that the Wigner transform is continuous 
in L 2 norm: 



E 



Ip2 ' 



< C / d£ sup 



o(t,v) V IE IIV'ilr euv'i - V2I 



(32) 

Thus the unitarity bound can also be used for the Wigner transform. The proof 
of the Boltzmann equation ED on an arbitrarily large kinetic timescale T uses an 
expansion up to order N ~ logt. The ladder terms give the gain term in the 
Boltzmann equation. The lowest order self-energy correction gives the loss term 
in the Boltzmann equation. It corresponds to the "gate" graph 



3.5 Long time scale: renormalization 

Because the ladder with n rungs is of order (X 2 t) n /nl, it diverges under diffusive 
scaling, and so do other graphs. To increase the time beyond A~ 2 , we need to do 
a renormalization. Formally, one can think of this as a resummation of the gate 
diagrams, which are of self-energy type, but this geometric series converges only 
for small A 2 i. A way to avoid such formal resummations is to change the way H 
is split into a "free" and an interaction part, i.e., expand around a different Hq. For 
£ > set 

6 e (a) = / dq ^— — . (33) 

J a - e(q) + is 

This is the value of the gate diagram at energy a in the Anderson model (in the 
QLM, the integrand contains an additional factor from the interaction function). 
The limit 0(a) = lim e ^o+ ©e(«) exists and is Holder continuous Q in a of order 
1/2. Let 

6(p) = e(e(p)). (34) 
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The idea is now to put \ 2 9{p) as a counterterm, which by construction subtracts ev- 
ery insertion of a gate diagram at the point a = e(p) where the particle propagator 
is singular. Because a and 7] appear only as auxiliary quantities in the expansion, 
it was necessary to take e — > above and define 6 in an a-independent way. 

The counterterm is added and subtracted so that the Hamiltonian does not 
change: let u(p) = e(p) + \ 2 9{p) and decompose 

H = uj{P) +U, U = XV - \ 2 6{P) (35) 

(where P denotes the momentum operator). The function u> can be thought of as 
a new dispersion relation of energy as a function of momentum. However, oj also 
has a negative imaginary part, roughly of order A 2 . More precisely, for d > 3 there 
is c > such that 

Im^p)<-cA 2 |H d - 2 . (36) 

Thus Hq is no longer selfadjoint. However, the negative sign of Im ui implies that 
the resolvent R a +i v is still well-defined, since the imaginary parts add up with 
the same sign. Correspondingly, the time evolution operator e~ lsH ° is no longer 
unitary but it remains bounded for s > 0. Both the Duhamel and the resolvent 
iteration are thus well-defined. Besides the new propagator (a + rq — o;(p)) _1 , the 
important change is that every factor U now also contains a counterterm insertion 
—X 2 6(p). The point about these iterations is that they can be stopped (or even mod- 
ified) after every expansion step. It is thus clear that one can group the counterterms 
that appear in the expansion together with the gates that get created when taking 
the average over the disorder. The cancellation among these two terms provides a 
small factor that makes such terms vanish in the diffusive scaling limit. Moreover, 
it is clear that one can implement rules for stopping the expansion independently 
of the subsequent disorder average. In particular, because the randomness is i.i.d., 
one can avoid moments beyond the power 2d by stopping the expansion when a 
given site has appeared in the collision history d times. The terms to which no 
such repetition or renormalization cancellation applies are expanded up to order 
\H\~ S ~ \- K ~ s , where 5 > depends on k. The intuition behind this is that 
certain graphs with n ~ \ 2 t ~ A~ K give the main contribution, and expanding up 
to an order that is A~ 5 higher leads again to small factors. 

The imaginary part of u) gives effectively a regularization 0(A 2 ) instead of 
0(?7) for the denominators, which changes the values of all diagrams significantly. 
In particular, the integral for one rung of the ladder becomes 

/ A 2 dp = 1 + Cq\ 1 ~°^ (37) 

(a-Lo(p+r)-irf) (/3-uj(p-r)+ir)) 
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where Co is a constant. Thus with this renormalization, the ladders become of 
order 1, so that one can go beyond kinetic scaling. Indeed, in the language of 
Feynman graphs, the main result can be stated informally as 

After renormalization, the sum of the ladder graphs for the Wigner 
transform converges to a solution of the heat equation in the diffusive 
scaling limit. 

The precise statements are Theorems 5.1, 5.2, and 5.3 in Ref.[7]. They involve 
in particular proving that the terms which do not correspond to pure up-down pair- 
ings vanish in the limit, and dealing with a number of technical complications 
which arise from the fact that one has to do an expansion to a finite order. 

3.6 The key estimate for controlling combinatorics 

We have had to leave out almost all technical details to avoid overloading the pre- 
sentation, but we should like to at least mention the heart of the proof here at the 
end, to clarify the main ideas about the Feynman graph expansion. 

Focusing on up-down pairings, we have to deal with a combinatorial problem 
of bounding the sum over the n! permutations a E S n . As mentioned, with an 
expansion to infinite order, one cannot get beyond the kinetic scaling because of 
this factor nl. The control of the remainders is done here by choosing an appro- 
priate stopping n for the expansion and by "beating down the combinatorics by 
power counting". That is, we prove exponential suppression of the values of Feyn- 
man graphs in the number of crossings they have, that is, loosely speaking, in their 
complexity. 

The precise notion capturing the complexity of a permutation a G S n is its 
degree d(a), defined as the number of non-ladder and non-antiladder indices. Es- 
sentially, the ladder indices are those for which a(i + 1) = + 1, and the 
antiladder indices are those for which a(i + 1) = a(i) — 1. 

Theorem 3.1 Let T a be the Feynman graph corresponding to a. There is 7 > 
such that for all a 

\Val(T a )\ <C\~< d{cT) . (38) 

This theorem is proven using a special integration algorithm for bounding the val- 
ues of large Feynman graphs [7 ]. 

The number of permutations with degree D is 

M n .D = \{a e S n : d(a) = D}\ < 2(2n) D . (39) 
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Expanding up to n = 0(X K 5 ), 5 > 0, we have by (l38l) . if 7 — k — 6 > 0, 

k k 

X^ d{a) = Y ^ d N n ,d < 2 (2X) d( ^ K ' S) < 0{X D ^- K - 5 ^). (40) 

^6S„ d=D d=D 

d(a)>D 

Thus the contribution from the sum of all terms with degree D > 2 is small if 
7 — k — 5 > 0, hence the essential restriction for the value of k is that of 7. As 
mentioned, one would hope to get close to 7 = 2 in (|38T >. but 7 has to be chosen 
smaller for technical reasons. 

4 Conclusion 

We have shown that, for random Schrodinger operators with a weak static disorder 
the quantum mechanical time evolution can be controlled on large space and time 
scales where a diffusion equation governs the behavior. The Schrodinger evolution 
is time-reversible - yet irreversibility on large scales emerges. This apparent con- 
troversy is resolved by noting that along the scaling limit microscopic degrees of 
freedom have been effectively integrated out. 

Although the expansion methods we use bear some resemblance to those of 
constructive quantum field theory, there are also a few noteworthy differences. 
First, because we analyze the time evolution at real time, the (near-)singularities 
of the propagators are located on hypersurfaces, and not at points, as would be the 
case in Euclidean field theories. The singularity structure is to some extent similar 
to that in real time Fermi surface problems, although there is no fixed Fermi surface 
here - the integrals over a and "test" all possible level sets of the function e{p), 
and this leads to a number of serious complications. Second, we are able to control 
the combinatorics of a straightforward Feynman graph expansion in momentum 
space, while the analysis in constructive field theory (to our knowledge, always) 
needs to be done by cluster expansions in position space to avoid the divergence of 
an infinite series of Feynman graphs. The reason for this is twofold: the unitarity 
bound allows us to do an expansion to a finite order, and our strong improvement 
(I38T ) over standard power counting bounds allows us to push this order so high that 
we can reach the scale where diffusion sets in, while still retaining control of the 
remainders. 

The genuine challenge is to show diffusion without taking scaling limits, i.e. 
for a fixed (small) disorder A and for any time independent of A. With expansion 
techniques, this would require to renormalize not only the self-energy to arbitrary 
order but also the four-point functions. Refining the self-energy renormalization 
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poses no fundamental difficulty. The correct renormalization of all four-point func- 
tions in this problem, however, remains a widely open problem. 
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